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i.  INTRODUCTION 


I 


A desirable  property  of  any  control  design  is  that  it  be 
insensitive  to  small  variations  in  the  parameters  of  the  controlled 
plant.  This  is  necessary  for  several  reasons.  For  instance,  the  mathe- 
matical model  can  only  approximate  the  physical  problem  so  that  the 
assumed  values  of  parameters  for  the  design  may  be  different  from  the 
actual  parameter  values  upon  implementation.  Also,  most  systems  suffer  from 
some  forms  of  unmeasurable  or  unpred ic tab le  variations  due  to  the  degen- 
eration of  physical  components  and  adverse  environmental  effects.  Hence, 
there  is  a need  for  general  design  methods  which  guarantee  some  degree 
of  insensitivity. 

In  this  thesis  the  reduction  of  sensitivity  of  feedback  systems 
to  parameter  variations  in  the  plant  is  investigated.  The  concept  of 
comparison  sensitivity  is  used  as  an  indication  of  the  sensitivity 
performance  of  a given  design.  In  particular,  the  question  of  sensitivity 
reduction  in  feedback  systems  which  use  state  observers  for  dynamic 
compensation  is  considered  leading  to  a design  procedure  which  guarantees 
sensitivity  reduction  with  respect  to  a particular  comparison  sensitivity 
reduction  criterion. 

The  problem  has  been  solved  for  the  full  state  feedback  case 
[1,2],  and  it  was  recently  solved  for  the  output  feedback  case  using 
dynamic  compensators  [3l.  In  both  instances  it  is  shown  that  sensitivitv 
reduction  is  directly  related  tc  some  optimal  control  law.  A short  over- 
view of  these  past  results  is  given  in  the  sequel  followed  by  a 


invariant  (LTI)  systems  which  are  state  controllable  and  state  observable. 


. COMPARISON  SENSITIVITY  WITH  FULL  STATE  FEED BACK 


l.l.  Definitions 

Comparison  sensitivity  relates  the  sensitivity  of  one  control 
scheme  to  the  sensitivity  of  another  nominally  equivalent  control  to 
identical  parameter  variations  in  the  controlled  plant.  The  potential 
benefits  of  using  state  feedback  to  improve  sensitivity  can  be  evaluated 
by  comparing  the  sensitivity  of  the  closed -loop  design  to  a nominally 
equivalent  open-loop  control.  The  development  of  these  concepts  has  lead 
to  the  definition  of  the  comparison  sensitivity  operator  [-»)  which  direct  Iv 
relates  the  open-loop  and  closed-loop  sens i t i vit ies , and  a sensitivity 
reduction  criterion  giving  sufficient  conditions  for  a particular  feedback 
control  law  to  guarantee  sensitivity  reduction  in  comparison  to  the  open- 
loop  control  [ 1] . 

To  make  the  above  discussion  precise,  consider  an  LTI  dynamic 
system  described  by  the  nth  order  state  equations 


xit)  • Ax  it)  t-  3u\t',  x i0)  “ 
xft)  £ Rn,  u it)  £ Rn* 


If  a state  feedback  law  is  given  by 


uit)  “ Fx  it ) . FtR" 


then  the  nominal  state  trajectory  is  given  by 
► \ . lA+BF't  o 

xn(t)  - e x . (3) 

Now,  consider  differentially  small  variations  in  the  parameters 
of  the  coefficient  matrices,  denoted  by  .'A,  5s,  and  5C.  Neglecting  second 


and  higher  order  terms,  the  variation  of  the  state  trajectory,  5x(t),  is 
given  as  the  solution  to 

5x(t)  **  5Ax (t)  + A5x(t)  + 5Bu (t)  + B6u (t ) , 5x(0)-0.  (4) 

Defining  the  open-loop  state  trajectory  variation  as 


5xo(t)  '"X^C'*‘(4)  with  5u(t)“0 
and  the  closed-loop  state  trajectory  variation  as 

:xc  (t)  * 5x(t)l(4)  with  5u(t)  ■ FSx  (t) 


(5) 


(6) 


leads  to  the  following  definition. 

Definition  (Comparison  Sensitivity).  Given  the  system  (1),  the  feedback 
law  (2)  defines  the  nominal  state  trajectory  (3).  With  (5)  and  (6) 
defining  the  open-loop  and  closed-loop  state  trajectory  variations,  the 
closed-loop  system  is  said  to  be  less  sensitive  than  the  nominally  equi- 
valent open-loop  system  with  respect  to  a weighting  matrix  Z,  if 


t 

b 


i 

5x^ (t)Z5xc  (t )dt 


5*o(OZSx  (t)dt 

0 


(7) 


for  all  t,  > 0.  The  matrix  Z in  (7)  is  a positive  semi -definite  symmetric 
nxn  matrix. 

The  comparison  sensitivity  matrix  is  a frequency  domain  operator 
defined  as  follows.  For  the  open-loop  and  closed-loop  systems  defined 
above,  it  can  be  shown  [2]  that  in  the  frequency  domain 


SXc(s)  = S(s)SXo(s) 


(8) 


where  5X  (s)  and  iX  (s)  are  the  Laplace  transforms  of  5x  (t)  and  5x  (,t) 

CO  CO 


Another  useful  form  for  Sys)  in  Is  given  by 


S(,si  “ Is-  ysI-A-BF^  8F. 


The  form  ^ 10 > is  derived  from  y 9 h by  applying  the  matrix  identitv 
(I+UV> " : * I-U(I+VT'>’IV. 


d.d.  State  Feedback  Results 

Having  defined  comparison  sensitivity,  the  following  results 
give  a direct  relationship  between  comparison  sensitivity  reduction  and 
an  optimal  control  which  minimises  a particular  quadratic  performance  i; 
Theorem  1 . [dl  A sufficient  condition  for  to  hold  for  all  r.  ^ d i< 

that  for  all  JtcX 

sVjortZSU*)  -2  £ 0 I 

where  3 (ja'i  “ S ys  V defined  in  yJi. 

s ■ j d? 

If  the  performance  index 

dO 

J " t . !!x  ' yt  iQx  yt > + u ' yt >Ru yt >'d t 

\ 

where  vj  * 2l0  and  R « R ' > 0 

is  defined,  it  is  well  known  that  the  state  feedback  law  which  minimises 
J is  given  by 

uyti  - -R  i3'\xyt''  v 

where  X is  the  unique  positive  definite  symmetric  solution  to  the  Riccat 


equation 


ka+a'k  + q-kbr~Vk  = 0. 

The  result  for  state  feedback  of  primary  interest  for  this 
research  may  now  be  stated. 

Theorem  2.  [ ll  For  a system  described  by  (1)  the  sensitivity  reduction 

criterion  (11)  is  satisfied  for  a feedback  law  (2)  if  there  exists  some  0 
and  R as  in  (12)  such  that 

F - -r'Vk  and  Z » f'bF.  (14 

In  other  words,  any  optimal  feedback  control  law  as  in  (13) 
reduces  sensitivity  in  the  sense  of  comparison  sensitivity  for  a weighting 
matrix  given  by  (14). 

It  should  be  noted  that  Z depends  integrally  upon  the  quadratic 
performance  index  which  is  optimized  and  is  therefore  not  at  the  direct 
disposal  of  the  designer.  Furthermore,  because  of  the  matrix  factor  R in 
(14),  the  rank  of  Z is  bounded  above  by  m,  the  number  of  inputs. 


«n  npiilL>w.ni  -J-‘  '-T'  ^ 


s 


J.  COMPARISON  SKNSiriVlVY  WITH  OCT  PIT  FKK1MIAOK 

As  In  other  areas  of  linear  system  theory,  assuming  a linear 
combination  of  the  states  rather  than  all  of  the  states  is  available  for 
measurement  greatly  complicates  the  feedback  control  sensitivity  problem. 

The  discussion  and  example  which  follow  reveal  some  of  the  reasons  why  the 
output  feedback  case  demands  theory  beyond  what  has  been  presented  foi  the 
state  feedback  case.  Some  recent  results  which  deal  with  this  more  genera’, 
situation  are  stated  and  these  results  form  the  basis  for  the  application 
of  state  observers,  which  is  the  contribution  of  this  research. 

3.1.  The  Need  for  Added  Dynamic s 

It  is  well  known  that  a feedback  law  which  is  strictly  a lineal 
t rans formation  of  the  outputs  of  a lineai  time  invariant  system  will  not, 
in  general,  be  sufficient  to  place  the  poles  of  the  system  arbitrarily. 
Hence,  reasonable  control  is  feasible  only  when  some  dynamic  compensation 
is  included  in  the  control  scheme.  The  added  dynamics  mav  be  m t he 
forward  path,  as  In  the  case  of  PID  control,  or  in  the  feedback  Loop,  as 
in  the  case  of  the  state  observer  implementation. 

Following  the  development  of  the  state  feedback  case  in  the 
previous  section,  it  would  be  reasonable  to  pose  the  comparison  seositivitx 
problem  in  terms  of  the  outputs  rather  than  the  states  and  to  determine  i 
sufficient  condition  undei  which  the  closed -loop  svstem  is  less  sensitive 
than  the  nominally  equivalent  open-loop  svstem.  The  following  example  shows 


that  even  for  a simple  case,  attempting  to  satisfy  such  an  output  sensitivitv 


TTT  X 


8 


criterion,  fails  for  conventional  design  procedures.  The  subsequent 
theory  presented  in  Section  3.3  establishes  that  any  added  dynamics 
become  inherently  involved  in  the  formulation  and  solution  of  the  sensi- 
tivity reduction  problem. 


3.2.  A Counterexample  Using  State  Observers 

Consider  the  following  conventional  design  procedure.  For  an 
unstable  plant,  an  optimal  state  feedback  control  is  found  using  the 
algebraic  Rtccati  equation.  By  Theorem  2,  such  a control  guarantees 
sensitivity  reduction  with  respect  to  a particular  weighting  matrix.  The 
state  feedback  is  implemented  by  direct  feedback  of  any  measurable  states 
and  the  remaining  unmeasurable  states  are  estimated  by  a full  order 
observer . 


As  an  attempt  to  extend  the  sensitivity  results  for  state  feed- 
back to  the  output  feedback  case,  an  output  sensitivity  reduction  criterion 
is  defined.  Using  the  above  design  procedure,  the  observer  dynamics  are 
designed  to  satisfy  the  sensitivity  criterion.  For  the  following  particular 
case,  such  a design  fails  to  result  in  a stable  observer. 

For  the  controllable  and  observable  LTI  system 

x(t)  ■ Ax(t)  + Bu(t),  x(t)  ■ x° 


y(t)  = Cx(t) 


U5) 


where 


0 

4. 

1 

* 

B - 

-1 

3 

l 

C * [0 


, i 


52 X' . 


the  state  feedback  control  which  minimizes 


9 


1 > r 2 

7 o -y  (c)  +u"  (c):dc 


(16) 


is  given  by 


u(t)  **  Fx  (t ) where  F =*  [2.88  -10.181 


(17) 


Following  the  approach  to  the  output  sensitivity  problem 

suggested  above,  suppose  the  sensitivity  reduction  criterion  to  be 

satisfied  is  given  by 
t 


j 6y“(t)dt  <u  5y~(t)dt 
o o 


(18) 


for  all  t^i-O.  The  second  state  will  be  fed  back  directly  and  the  first 
state  will  be  estimated  using  a full  order  state  observer  of  the  form 

x = [a  + LC] x + Bu  - Ly  (19) 

where  A,  B,  and  C are  defined  by  the  system  (15)  and  L is  the  observer 
gain, 


l. 


(20) 


which  is  chosen  to  satisfy  (13)  and  to  stabilize  (19). 

It  can  be  shown  that  a necessary  and  sufficient  condition  for 
(13)  to  hold  for  this  design  is  (see  Appendix) 

s'*  - (8 .88+i0 )s3+  (28.04-1.88/^+5  .83/.,  )s“-  (25  .-+4-5 .64/ ^+6 .64/.,  )s 
-24. 8+3. 76/, -5.76/, 


1 


s *(1.30-Z2)s  -(12.68-/  +4. 30/, Js1- (15. 26-4. 30/ j +2.22/.,) s 

-4.44+2.22/, 


<1  (21) 


1 


where  s = jit,  for  all  i>  € R 


which  is  true  if  and  only  if 


. • ' . ; 


■■host. 


10 


r..5 


d -(28.04-1.882L+5.88^o)a)“-24.8+3.76iL-5.76.i9]“+  (8.38+2^ 
-(25.44-5.6421+6.642,)a']2-[a;i++(12.68-xJ4.302.,)uu2-4.44+2.22ii:2 
-[  (1.30-22)oj3  + (15 .26-4.302L+2.222o)a']2  < 0. 


(22) 


Expanding  these  terms  and  collecting  the  coefficients  of  the  powers  of  d 
leads  to  the  highest  power  of  <jq  being  six  with  the  coefficient 


-4.28  + 5.762^. 


Hence,  a necessary  condition  for  (22)  to  hold  is 


21  < .74. 


However,  stability  of  the  observer  requires  that 


l >2. 

1 im 


(23) 


(24) 


(25) 


Therefore,  the  conditions  (24),  and  (25)  imply  that  no  observer  imple- 
mentation will  satisfy  both  the  sensitivity  criterion  and  the  stabilization 
requirement . 

This  simple  example  shows  that  added  dynamics  lead  to  sensitivity 
problems  if  the  output  sensitivity  criterion  (18)  is  to  be  satisfied.  In 
the  following  section  a general  theory  for  comparison  sensitivity  with 
output  feedback  will  be  presented  which  shows  that  the  added  dynamics  of 
a compensator  must  be  included  in  any  comparison  sensitivity  reduction 
criterion. 

3.3.  Output  Feedback  Results 

This  section  is  based  entirely  upon  the  theory  developed  in  [3]  . 

To  study  comparison  sensitivity  for  the  output  feedback  case 
the  sensitivities  of  all  variables  involved  in  the  feedback  law  must  be 


Lt 


11 


included  in  Che  sensitivity  criterion.  Hence,  when  states  from  a dynamic 
compensator  as  well  as  the  process  outputs  are  fed  back,  both  the  plant 


outputs  and  the  compensator  states  must  be  used  in  the  comparison  sensi- 
tivity analysis. 

To  make  this  explicit,  consider  the  LTI  controllable  and  observ- 
able system  given  by 

x(t)  “ Ax (t ) + Bu  (t ) , x(0)  = x° 

(26: 

y(t)  » Cx  (t) 

x(t)cRn,  y(t)  (r‘,  u(t)=Rm. 

Assume  rank(B)*=m  and  rank  (C)  =1  with  m<n  and  2<n.  Suppose  the  control 
u(t)  is  given  by 

u(t)  - Hz(t) +Ky(t),  H € RIWk , K € r"'*2  (27] 

where  z(t)  is  the  state  of  the  kth  order  dynamics  of  the  compensator  given 


z(t)  = Fz(t)  + Gy  ( t ) + Du ( t ) , z(0)  = z°,  zcR^.  (2i 

The  constant  matrices  F,  G,  and  D are  of  appropriate  dimensions. 

The  system  (26),  compensator  (23),  and  the  control  (27),  define 
a composite  system  with  states  x(t)  and  state  equations 


x(t)  ■ A x(t)  + Bu(t)  x(0) 

y(t)  - cx(t) 

u(t)  'iCy(t) 


u_ 


mv  i,mm.rr 


' : '' r ’ \ . 
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where 


F 

GC 

D 

A = 

B = 

_0 

A . 

B 

I. 

0 

c - 

k 

K “ 

Ch 

0 

C _ 

For 

differential  variations 

and 

5C, 

the 

resulting  first 

output  variables  satisfy 

5x(t)  = 5a  x (t)  + A$x(t)  + 5Bu(t)  + B6u  (t ) , 5x(0)  = 0 . 


(30) 


In  a manner  analogous  to  the  state  feedback  case,  the  open-loop  output 
trajectory  variation  is  defined  as 


5y0(t)  = 5y(t)|(30)  with  5u(c)  = 0 


and  the  closed-loop  output  trajectory  variation  is  defined  as 

sy„(t)  - 5y(t)|  „4i.u 


(31) 


(30)  with  0u (t)  = K5y (t) ' 


(32) 


The  comparison  sensitivity  criterion  [4]  may  now  be  defined. 
Given  the  composite  system  (29)  with  input  u(t)=K7"(t)  defining  the 
nominal  output  trajectory,  the  closed-loop  system  is  less  sensitive  than 
the  nominally  equivalent  open-loop  system  with  respect  to  a positive  semi- 
definite  symmetric  (n+k)X(n+k)  weighting  matrix  Z,  if  for  all  t^>0 


fcl  *1 

J 5y,!(t)Z5y  (t)dt  < „ 3y^(t)z5yo(t)dt. 
o o 


(33) 


o 


i 


1 


In  a manner  analogous  Co  the  state  feedback  case  [. 4- , 2 ] , a 
sufficient  condition  for  (33)  to  hold  is  given  in  Theorem  3. 

Theorem  3.  [4]  The  plant  and  compensator  defining  the  composite  system 

(29)  satisfy  the  condition  for  comparison  sensitivity  reduction  (33)  if  for 
all  m in  R 

[l4C(-jtuI-A-BK  C)'1  B Kj  'zCl+C(ju.'I-A-3  K C)"1  B KJ-Z  < 0.  (34) 

Letting  (34)  define  the  comparison  sensitivity  reduction 
criterion,  the  following  theorem  parallels  the  result  of  Theorem  2 for  the 
state  feedback  case  by  relating  the  satisfaction  of  (34)  to  the  solution 
of  an  optimal  control  problem. 

Theorem  4.  [3]  For  a plant  and  compensator  given  by  the  composite  system 

(29)  with  m<k+A,  a particular  parameterization  of  the  compensator  (23) 
and  (27)  satisfies  the  sensitivity  reduction  criterion  (34)  for  a positive 
semi-definite  symmetric  weighting  matrix  Z with  rank(Z)=m  if  and  only  if 
the  following  three  conditions  are  satisfied 


rank(K)  = m 

(35) 

_ 

Z = K RK 

(36) 

-1  — 1 

KC  =*  -R  BP 

(37) 

where  P is  the  positive  definite  symmetric  solution  of  the  n-Hc  order  Riccati 
equation 

PA  + A ' P + Q-  PBR*  L3 ' P - 0 

for  some  positive  definite  symmetric  matrices  Q and  R of  appropriate 
dimensions . 


r 
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The  conditions  given  by  (36)  and  (37)  may  be  viewed  as  a direct 
extension  of  the  state  feedback  result  since  K C is  actually  the  state 
feedback  matrix  for  the  composite  system  (29),  and  (37)  demands  that  it 
satisfy  an  optimal  control  law. 

In  [3],  the  conditions  of  Theorem  4 are  applied  to  a canonical 
form  for  the  composite  system  which  leads  to  a design  procedure  for 
constructing  a dynamic  compensator  which  satisfies  the  sensitivity 
reduction  criterion.  The  details  of  these  results  will  not  be  discussed 
here,  but  those  points  which  are  used  in  the  development  of  a design 
procedure  using  observers  will  be  presented  briefly. 


Lemma  1 . (Plant  Canonical  Form).  Given  a plant  described  by  (26)  with 
rank(CB)=m^,  there  exists  nonsingular  matrices  T,  nXn,  and  U,  mXm,  such 


CT  = [ I, 


where 


T~ 1 B U 1 


B*  - 
> D"> 


1 €Rtn'^Xni. 


This  lemma  clearly  follows  by  construction  and  implies  that 
there  are  bases  for  the  state  and  input  spaces  such  that  (38)  and  (39) 


hold. 


It  is  shown  in  [ 3 1 that  a necessary  condition  for  rank(Z)=m  is 
that  rank(,C3)*  m.  Bv  selecting  D in  the  compensator  to  satisfy  this  rank 


1 


condition  it  can  be  shown  that  there  exist  nonsingular  matrices  T_, 

(k+n) X (k+n),  and  T , (k+Z)X (k+4 ) , such  that  the  composite  svstem  can  be 
o 

transformed  to  what  will  be  called  the  composite  system  canonical  form 


given  by 


A A - -1 
A - T AT 
s s 

„ . r i 

„ a „ — I m 

B = T B = l n 

s L 0 . 


C = T ct'1  = [ I , . 0]  . 

o s k+Z 


With  (A,B,C)  defined  above,  the  following  equivalent  conditions 
for  sensitivity  reduction  can  be  established. 

Theorem  5.  [3]  For  the  composite  system  (29)  with  the  plant  matrices 
in  plant  canonical  form  the  conditions  (35),  (36),  and  (37)  of  Theorem 
are  equivalent  to 


rank  (C  B)  = m 


and  the  eigenvalues  of 


? ~ r _ ,,  - „mX  (k+i-m) 

A22  +A21W“Ik+Z-m  ’ W tR 


have  negative  real  parts,  where  A is  any  matrix  similar  to  A defined  in 
(40)  and  A is  partitioned  as 


t . Al1  A12  ; .-j*  I t o (n*Ht-m)X  (n+k-m) 

A a » , A,.  »K  , 

A A 
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The  proof  of  this  theorem  given  in  [3]  is  constructive  in  that 
it  leads  to  the  following  design  procedure.  A dynamic  order  k is 
selected  for  the  compensator  and  the  compensator  is  designed  with  D,  F, 
and  G as  in  (28)  and  satisfying  (41).  The  composite  system  is  trans- 
formed to  the  composite  system  canonical  form  and  a matrix  W is  sought  to 
satisfy  condition  (42)  of  Theorem  5.  If  no  such  matrix  exists,  the  order 
k is  increased  until  such  a W can  be  found.  It  is  shown  in  [3]  that  for 
some  k<n-2+ra-m,  there  can  be  constructed  a compensator  and  W to  satisfy 
(42). 

Once  W has  been  found,  the  following  matrices  are  defined 

Oil 


w =* 

I 

m 

W!  I,  .a 
K-nc  -m 

0 

^k+n-m 

^ — 

A = 

A11 

A12 

- W 

1 

> 

ro 

»— * 

A22 

(44) 


(45) 


A positive  definite  symmetric  matrix  Q,  (k+n-ra)X  (k+n-ml , is  selected  and 
the  Lyapunov  equation 


PA^  + A22P  + Q - 0 

is  solved  for  P.  Since,  by  (45)  .A,.,  is  equal  to  che  matrix  (42),  the 
matrix  P is  a positive  definite  matrix  uniquely  solving  (-+6).  Now,  an 
mXm  positive  definite  matrix  R is  selected  such  that 


(46) 


4 -il2slZ'u  > o 


(47) 


and  a (k+n-m)\ (,k+n-m)  positive  semi -definite  symmetric  matrix  V is  chosen 
such  that 


Q*  ^ V + a!lPq'1PA^1  > 0. 


The  Riccati  equation 


where 


P*A*  + A*P*  - P*R*P*  + Q* 


A Au+A12Q  PA21 


is  then  solved  for  the  unique  positive  definite  symmetric  matrix  P* 

(such  a solution  is  guaranteed  by  the  selection  of  V and  R) . The  composite 
system  is  then  stabilized  and  the  sensitivity  criterion  (34)  is  satisfied 


if  K is  given  as 


-1  *r 

K - -R  P*L I -W^T  . 

o 


As  was  mentioned  above  the  details  and  analysis  of  this  design 
procedure  are  given  in  [3].  The  problem  of  selecting  the  dynamic  order 
of  the  compensator,  k,  and  the  matrices  D,  F,  and  G is  solved  in  the  next 
section  by  implementing  state  observers  as  the  compensator  dynamics.  This 
leads  to  a comprehensive  design  procedure  and  some  simplification  of  the 
interpretation  of  the  algorithm  through  the  selection  of  the  transforma- 
tion matrices  T and  T . 

o s 


.. . 
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•+.  COMPARISON  SENSITIVITY  WITH  STATE  OBSERVERS 

4.1.  Whv  State  Observers? 

The  results  presented  in  the  last  section  are  a solution  to 
the  problem  of  comparison  sensitivity  reduction  when  output  feedback  is 
used.  It  is  not  a totally  satisfactory  solution,  however,  because 
although  the  sensitivity  criterion  is  satisfied  and  the  system  is 
stabilized  by  the  compensator,  there  is  no  explicit  form  for  the  dynamics 
of  the  compensator  and  there  is  no  direct  control  over  the  eigenvalues  of 
the  overall  system. 

It  is  reasonable  to  ask  how  known  dynamic  compensation  techniques 
fit  into  the  sensitivity  design  specifications.  Hence,  the  purpose  of  this 
research  was  to  study  the  sensitivity  characteristics  of  feedback  systems 
which  use  state  observers  for  the  compensator  dynamics.  The  use  of  state 
observers  is  desirable  since  their  structure  is  well  known  and  the  design 
of  observers  is  easily  computerized.  More  importantly,  it  is  a property 

of  systems  which  have  state  observers  that  the  observer  poles  are  also 

poles  of  the  overall  feedback  system  and  determine  the  speed  with  which 
the  error  of  the  state  estimation  diminishes.  Finally,  since  many  feed- 
back laws  are  given  in  terms  of  state  feedback  and  the  states  often  have 

physical  significance,  it  is  desirable  to  have  available  an  estimate  of 
the  unmeasurable  states  to  implement  a state  feedback  scheme  and  to  monitor 
the  performance  of  the  system. 


19 


TjT 


4.2.  The  Full  Order  Observer 

The  compensator  dynamics  (28),  if  chosen  to  be  a 
observer  (i.e.  k = n),  will  have  the  form 


full  order 


F =*  A+LC 

G =*  -L  (51) 

D ■=  3 


where  the  nxi  matrix  L is  chosen  such  that  the  eigenvalues  of  F will  have 
negative  real  parts.  Usually  the  eigenvalues  of  the  observer  are  chosen 
to  be  much  faster  than  the  remaining  feedback  system  eigenvalues  so  that  the 
estimation  error  goes  to  zero  relatively  quickly. 

Assume  the  system  (26)  is  in  the  plant  canonical  form  given  in 
Lemma  1.  It  is  clear  that 


rank  (C  B) 


rank 


(52) 


so  that  the  rank  condition  (41)  of  Theorem  5 is  satisfied. 

Now  partition  the  matrix  A as  in  Figure  1 and  define  T and  T 

o s 

as  in  Figure  2.  It  is  simple  to  verify  that  the  nonsingular  matrices  Tq 
and  Ts  will  transform  the  composite  system  (29)  for  the  plant  (26)  and 
observer  defined  in  (51)  to  composite  system  canonical  form.  The  explicit 
form  for  A (40)  in  this  case  is  given  in  Figure  3.  With  A partitioned  as 
in  (43),  the  stability  of  the  matrix  (42)  is  shown  in  Figure  4 to  depend 
upon  the  stability  of  the  observer  and  the  existence  of  a matrix  W which 


siabilizes  the  pair 


r 


All  • ‘ ' A14 


aii€r 


X 


A33€r 

v£r 


(m-nij)  X (m-m^) 


(n-£+m^-m)  X (n-f+m^-m) 


Figure  l.  Partition  of  Plant  Matrix,  A. 


0 0 


I 0 

m-ra^ 


n-/  -ra-Hn. 


I 0 


! ' ~Z] 


(n-i)K  { In-i  °(n-/)xi 


Figure  2.  and  Tg  for  Full  Order  Observer  Case. 


2 3 


(53) 


The  existence  of  such  a mateix  W is  guaranteed  by  Che  following  lemma. 
Lemma  2 ■ Given  the  controllable  pair  (A,B)  where  B is  of  the  form  given 
in  (39),  if  A is  partitioned  as  in  Figure  1,  then  the  pair  (53)  is 
controllable . 

Proo f . Let  \ ^ , i " 1,2, . . . ,a  be  the  n eigenvalues  of  A.  The  pair  (A,B) 
is  controllable  if  and  only  if  [5] 


rankrA-X^I;  B]  = n for  i = l,2,...,n. 
Since  B is  of  the  form  in  (39),  (54)  implies 


An"XiI 

A12 

A13 

A14 

1 

1 

I 

ml 

0 

A21 

A22AiI 

A0  ^ 

A24 

1 

| 

0 

0 

A31 

A32 

A33*X i1 

A34 

1 

1 

0 

I 

m-m. 

A,  , 

L * L 

A,  0 
^2 

A43 

A, . -\  .1 
•44  i 

1 

0 

0 

(54) 


(54a) 


which  is  true  if  and  only  if 


rank 

A21 

A22_Xi*  A23  A24 

=*  n-m. 

(55) 

_A4 1 

A / /%  A/<n  A/  , * V i X 

^♦2  4J  *4*4  I J 

Rearrangement  of  the  matrices  in  (55)  gives  the  contro 1 lab i 1 i ty  of  the 
pair  (53). 

From  Figure  4 and  Lemma  2 it  is  clear  that  there  exists  a matrix 
W such  that  condition  (42)  of  Theorem  5 is  satisfied  for  any  stable  full 


order  observer. 


■—  ~ « — r ' — z 
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Therefore,  Che  following  design  procedure  will  lead  to  a dynamic 
compensator  which  will  guarantee  comparison  sensitivity  reduction  with 
respect  to  some  weighting  matrix  Z of  rank  m.  Select  a full  order  state 
observer  (51)  and  choose  W so  that  the  pair  (53)  is  stabilized,  then 
follow  the  design  algorithm  given  in  the  previous  section.  The  design 
algorithm  will  be  given  explicitly  in  Section  5.1. 

One  important  property  of  dynamically  compensated  systems  which 
have  full  order  state  observers  in  the  feedback  loop  is  that  the  eigen- 
values of  the  observer  are  also  eigenvalues  of  the  overall  feedback  system. 
This  can  be  demonstrated  by  applying  the  similarity  transformation 


T 


-I 

I 


(56) 


to  the  composite  system  which  has  full  order  observer  dynamics 

fA+LC+BH  -LC+BH 1 


x(t)  - A^x(t)  where  Af 


BK 


A+BK 


for  arbitrary  output  feedback 

u(t)  = K y (c)  = Hz  (t ) + KCx(t). 


The 


similarity 

TAfT~ 


transformation  (56)  gives 
A+LC  0 

a 

BH  A+B (KC+H) 


(57) 


(53) 


(59) 


Note  also  that  the  other  n eigenvalues  of  the  system  depend  only  upon  the 
selection  of  K,  the  output  feedback  law. 


TITT 


rrr  x . 


■MW".  __ 


— t- 
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4.3.  The  Reduced  Order  Observer 

Again  assume  that  the  plant  matrices  (A,B,C)  are  in  plant 
canonical  form.  The  compensator  dynamics  (28)  if  chosen  to  be  a reduced 
order  observer  (i.e.,  k = n-2)  have  the  form 


F - Aj2  + Lk*u 
G - -FL  + A*l  + LA*l 
D = B*  + LB* 


A11 

A12  l 

A13 

A14 

A*  1 A*  1 

Afi  | A12 

B 

A21 

^2  | 

A23 

*24 

a*  A* 

21  1 22 

A31 

' a"  r A"  ” 
32  I a33 

A34 

_A41 

A42  1 

A43 

A44_ 

n* 


are  defined 

in 

(39) 

and 

_ 

(m-m^)Xm1 

hi 

li2" 

L11 

€ R 

L * 

> 

(k+m. ~m)X (2 -m? ) 

L_. 

T.,  , 

T._  . 

<=  R 

L ^ 

22J 

22 

• 

(60) 


L is  chosen  such  that  F has  eigenvalues  with  negative  real  parts. 

Clearly,  with  the  dynamics  of  the  compensator  defined  by  (60) 

we  have 


rank (C  B)  a rank 


D 

CB 


m 


(61) 


so  (41)  in  Theorem  5 is  satisfied. 

Now  partition  A as  in  Figure  1 and  define  T and  T as  in 

o s 

Figure  5.  It  can  be  verified  that  T and  T as  given  will  transform 

o s 

the  composite  system  to  the  composite  system  canonical  form.  A is  given 





in  Figure  6 and  Figure  7 shows  that  the  stabilization  criterion  (42)  in 
Theorem  5 is  satisfied  for  any  stable  reduced  order  observer  and  W chosen 
to  stabilize  the  pair  (53).  Hence,  by  Lemma  2,  a matrix  W exists  as 
required  for  any  choice  of  the  observer  dynamics. 

The  design  procedure  using  reduced  order  observers  is  identical 
to  that  using  full  order  observers  except  that  W is  entirely  determined 
by  stabilizing  the  pair  (53).  Also,  the  property  of  the  observer  eigen- 
values being  eigenvalues  of  the  composite  feedback  system  holds.  This  is 
shown  bv  writing  the  composite  feedback  as 


where 


x(t)  = Af x ( t ) , x(t) 


z(t) 

xL(t)  , 


xL(t)  cr" 


F+DH 

G+DH 

0 

3*H 

Ajl+BfH 

A* 

12 

3*H 

A* 

(62) 


and  applying  the  similarity  trans formation 


(63) 


which  gives 


TAfT_1 


F 0 

L BH  a+B'Ik+hl:h.J 


(6a) 


From  (64)  it  is  noted  that,  unlike  the  full  order  observer  case,  the 
remaining  n eigenvalues  of  the  feedback  system  depend  explicitly  upon  the 


observer  matrix  L. 


I 


. . . . Some  Co r o 1 1 n r los 

The  results  of  the  preceding  two  sections  can  be  summarized  as  a 
corollary  to  Theorem  ^ . 

Corol '. ary  l.  Given  an  LTI  system  described  by  with  coefficient 

matrices  (A.B.O  in  plant  canonic  form,  for  arbitrary  state  observer 
dynamics  of  full  order  (.51)  or  reduced  order  (.60 ) there  exist  output  feed- 
back gains,  K,  such  that  the  composite  system  - g ^ is  stable  and  insensi- 
tive to  parameter  variations  in  the  plant  with  respect  to  a comparison 
sensitivity  weighting  matrix  Z.  given  by  with  mti  ■ Kv^t). 

Furthermore.  the  eigenvalues  of  the  observer  are  eigenvalues  of  the 
composite  feedback  system. 

In  ( J ] it  is  shown  that  there  exists  a dynamic  compensator  with 
order  k<n-i-Hn-m,  such  that  the  sensitivity  reduction  criterion  is 
satisfied.  The  results  for  the  reduced  order  observer  give  the  following 
improvement  on  the  upper  bound  of  the  dynamic  order  of  the  compensator 
needed  to  satisfy  the  sensitivity  reduction  criterion. 

Corol'.  ary  1.  There  exists  a dynamic  compensator  which  satisfies  the 
conditions  of  Theorem  b for  some  dynamic  order  k < n - < . 

Proo  f . It  has  been  shown  that  the  conditions  of  Theorem  5 can  be 
satisfied  using  a reduced  order  observer  as  the  dynamics  for  which  k ■ n-t  . 

The  particular  choices  of  the  trans format  ion  T and  T for  both 

O S 

full  and  reduc.- ' order  observer  cases  lead  to  a useful  mterpretat  ion  of 
the  sensitivity  weighting  matrix.  Writing  out  v In N using  the  form  of  k 
given  in  (50)  gives 


w 


where 


Z = l'  HT 

0 o 


p*r"  ~p* 


-p*r"1p*w 


w'p*r"lp*w 


Consider  the  full  order  observer  case  with  as  in  Figure  2.  If  the 
output  y is  decomposed  as 


v * Lzj  z.[  z!  y ^ y^j 

m-m^ 

z^  £ R , z,  € R , z^  € R 

n-i+m, -m  m^  £-m, 

z.  £ R L , y.  £ R , v0  £ R 
4 I "2 


V " [y[  z3  2 4 y2!  2{-yl:  z2_y2" 


Since  the  states  of  the  observer  are  estimates  of  the  states 

of  the  plant  and  y(t)  is  the  first  l states  of  the  plant  when  C is  given 

by  (38),  the  last  two  components  of  the  vector  in  (67)  are  the  error  of 

the  estimation  of  the  state.  Also,  the  second  and  third  components  are  the 

estimates  of  the  n-i  unknown  states. 

Similarlv  for  the  reduced  order  observer  case  with  T as  in 

o 

Figure  5,  and  y decomposed  as 

m-m^  n-i+m^-m 

y ‘ -zl  z2  yl  y2-  ’ > z,  “ R , 

(68) 

m1  .c-m. 

v ^ £ R , v^£r 


V * Lyl-  zl'LllyrL12y2  : z2'L21yrL22y2  : y2 ' ' 
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The  design  of  the  reduced  order  observer  gives  the  estimate  of  the  n-Z 
unmeasurable  states  by  z-Ly  where  y is  the  Z measurable  states.  Hence, 

(69)  shows  that  T^y  gives  the  measurable  states  and  the  estimation  of  the 
unmeasurable  states. 

Since  the  states  of  a system  often  have  physical  significance, 
the  following  corollary  which  summarizes  the  above  discussion  may  be  useful 
in  evaluating  the  sensitivity  reduction  capabilities  of  a particular 
weighting  matrix  resulting  from  (36). 

Corollary  3.  Design  of  the  feedback  gains  according  to  the  state  observer 
schemes  defined  by  Corollary  1 leads  to  a sensitivity  weighting  matrix  H 
(65)  which  weights  the  known  states,  estimates  of  the  unknown  states  and 
in  the  full  order  observer  case,  the  error  of  the  state  estimation  given 
by  (67)  and  (69). 


I 


i 
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5.  THE  DESIGN’  ALGORITHM 

A summary  of  the  design  procedure  described  in  the  discussion 
following  Theorem  5 and  extended  by  the  use  of  state  observers  in  the 
preceding  section  will  be.  presented  in  the  following.  An  interactive 
computer  program  was  written  to  implement  the  design  algorithm  and  to 
enable  adjustment  of  the  arbitrary  parameters  iteratively  to  obtain  the 
desired  overall  design  objectives.  This  program  is  discussed  briefly  in 
conjunction  with  a descriptive  flowchart.  An  example  of  the  steps  involved 
and  comparison  to  a conventional  design  technique  is  given. 

5.1.  Flowchart  and  Interactive  Software  Description 

Figure  8 ?s  a flowchart  of  the  design  procedure  and  shows  the 
points  at  which  redesign  of  parameters  can  take  place.  Each  step  of  the 
algorithm  will  be  discussed  briefly  with  the  numbers  referring  to  blocks 
on  the  flowchart. 

1.  The  basic  system  parameters  must  be  determined  and  the  plant 
canonical  form  is  assumed.  Hence,  before  entering  a general  plant  descrip- 
tion, the  triple  (A,B,C)  must  be  transformed  to  the  plant  canonical  form. 
The  computer  accepts  the  A matrix  and  constructs  the  B and  C matrices 
according  to  the  dimension  specifications 
n - dynamic  order  of  the  system 
m - number  of  inputs 
l - number  of  outputs 
m,  - rank  of  CB. 
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Begin 


LYAP 


Start 


Accept  Plant  j 

2 \Cont80bs?V— 


Yes 


Place  n-m  Poles  with  F* 


4<I°lE>^ — 

[F 

5 Place  N Poles.T0,Ts 
6 1 W=[F*W2] 

I 


° A = W'1  Ts  ATs-1  W 

9 P from  Lyapunov  Eqn. 

10  p*  from  Riccati  Eqn. 

I 

11  Feedback  System  Calc. 

12  Display  Final  Results 

13  ^Continue  ?^> — 

1 Yes 

f Go  to  Selected  A 
V Re-entry  J 


Terminate 


■ Riccati  Design 
T0,Ts,W  = F* 


-(  S|°P  ) 


f P-5878 


Figure  8.  Design  algorithm  flowchart. 


2.  The  control labi lity  of  the  pair  (A,B)  and  the  observability 


of  the  pair  (A,C)  are  tested  and  the  program  terminates  if  either  test 


fai  Is . 


3.  The  stabilization  of  the  pair  (53)  is  accomplished  using 


pole  placement  of  the  n-m  poles  available.  The  resulting  feedback  matrix  F* 


is  used  to  determine  W after  the  observer  order  has  been  selected. 


4.  Reduced  or  full  order  observer  design  is  selected. 


5.  For  the  full  order  observer  the  matrix  L is  determined  using 


pole  placement  to  stabilize  the  n eigenvalues  of  F in  (51).  T and  T are 

o s 


constructed  as  specified  in  Figure  2. 


6.  The  matrix  W2  is  arbitrarily  chosen  by  the  designer  and 


V - IF*  W2]  as  specified  in  Figure  4. 


7.  For  the  reduced  order  observer  the  n-2  order  dynamics  are 


stabilized  using  a Riccati  optimization  design  which  requires  a 


(,n-£)X(n-£)  positive  definite  matrix  Q and  an  Zx.l  positive  definite  matrix 


R.  L is  then  found  as 


L - -KA^R*1 


where  A'p  is  defined  in  (60)  and  K is  the  (n-f)X(n-i)  unique  positive 


definite  symmetric  solution  to  the  algebraic  Riccati  equation 


. ,*  .-  , „ . v\*'d"  V,* 


KA;.,  + A",K  + Q - KApR  ApK  = 0 


Ts  and  T ^ are  constructed  as  in  Figure  5 and  W * F*  as  specified  bv 


Figure  7 and  Lemma  2. 


I 


3.  Set  up  A as  in  (45). 
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l> . Acco.pt  Q and  solve  for  1’  in  the  ksn-m  dimensional  Lyapunov 
equation  (4b'). 

10.  Accept  R ^ nut  V as  in  (-*7)  and  (•sS’'  and  solve  the  m\m 
Riceati  equation  (49)  for  P* . 

11.  The  feedback  gain  K (501  can  now  be  calculated  as  well  as 
the  sensitivity  weighting  matrices  Z (dbl  and  H (t>51.  The  resulting  feed- 
back system  [(,571  or  (t'Jl|  is  constructed  and  its  eigenvalues  are  calculated. 

17.  All  pertinent  results  of  the  final  solution  are  displaced. 

13.  The  program  can  he  terminated  or  re-entered  at  points 
indicated  on  the  flowchart  to  redesign  some  ispect  ot  the  system. 

In  summary,  the  design  procedure  involves  the  solution  of  the 
observer  feedback  gains  (nth  order  or  (n-ilth  order'l,  the  solution  of  a 
vn+k-  mi  order  Lyapunov  equation  and  an  mth  order  Rteeati  equation. 

Obviously,  there  are  myriads  of  free  parameters  which  are  to  he  selected 
bv  the  designer  and  it  may  take  several  design  attempts  to  accomplish 
satis  far  t o ry  re suits. 

Many  of  the  mat hemat ica l computer  subroutines  used  were  taken 
from  the  LINSYS  package  In]  which  was  also  used  to  make  the  plots  for  the 
fo  l low ing  example . 


3.3.  An  F.\anvp_le 

As  an  example  of  the  design  procedure,  consider  the  following 
model  of  the  longitudinal  dynamics  of  an  aircraft  [7],  In  Table  l the 
seven  states  md  three  inputs  are  defined.  Tlte  first  five  states  come 
from  the  simplified  model  of  the  dynamics. 


The  additional  two  states  are 


defined  as  integral  controls  to  eliminate  steady  state  error  for  constant 
set  points  of  velocity  and  pitch  angle.  Following  the  design  in  [7],  it  is 
assumed  that  the  velocity,  pitch  angle  and  altitude  are  measured  directly 
and  the  pitch  angle  rate  as  well  as  the  angle  of  attack  must  be  estimated. 
In  other  words,  an  observer  must  be  used  to  estimate  the  states  one  and 
four . 


Using  the  linearized  coefficients  and  flight  conditions  given 
in  [7]  the  state  equation  in  the  form  (26)  has  the  coefficient  matrices 
(A,B,C)  given  in  Figure  9.  In  Table  2 certain  parameters  defining  the 
flight  conditions  are  given  along  with  perturbed  values.  Assuming  these 
perturbations,  the  resulting  A and  B matrices  are  given  in  Figure  10. 

The  C matrix  is  unperturbed  since  the  states  are  measured  directly. 

To  implement  the  design  procedure,  the  plant  matrices  (A.B.C) 
were  transformed  to  the  plant  canonical  form  as  in  Lemma  1 by  the  non- 
singular matrices  T and  U in  Figure  11.  The  transformed  nominal  A matrix 
is  shown  in  Figure  12.  Note  that  under  the  similarity  trans formation  T, 
the  state  variables  defined  in  Table  1 are  arranged  as 


Lx2  x3  x5  X6  x7  X1  x4- 


(-2) 


The  design  algorithm  was  used  iteratively  with  the  final  choices 
of  free  parameters  and  resulting  observer  and  feedback  tutrices  given  in 
Figure  13. 


For  a comparison  of  the  sensitivity  reduction  properties  of  the 
design  algorithm  to  a conventional  design  procedure  usi.ng  a reduced  order 
observer,  the  feedback  gains  were  calculated  assuming  the  state  feedback 


mmm 


Table  l.  Definition  of  State  Variables  and  Inputs 

For  Aircraft  Control  of  Longitudinal  Dynamics 


X1 

a 

angle  of  attack 

X6 

" x0  - [Reference  Velocity  ~ 100] 

*2 

a 

velocity  a.  100 

x7 

* - [Reference  Pitch  Angle] 

X3 

m 

pitch  angle 

U1 

• elevator  angle 

x4 

a 

pitch  angle  rate 

U2 

= flap  angle 

X5 

at 

altitude  100 

U3 

* throttle 

Table  2.  Flight 

Condition  Parameters 

Variable 

Nominal  Value 

",  of  Perturbation 

Perturbed  Value 

Center  of 
Graviey 

.20  ft 

+ 5" 

.21  ft 

Velocity 

190.66  ft/sec 

+ 5" 

200.19  ft/sec 

We i gh  t 

4000  lbs 

- 5% 

3800  lbs 

Air  Density 

.004342  slugs /ft 3 

+ 5" 

.0050S4  slugs/f 

Moment  of 
Inertia 

2050  slugs-ft“ 

**  «-  a 

'■> 

2009  slugs-ft“ 

fc---  -'i- ... 


■ m . 


■ ■ ' . ' - mmm  mm  bm — ■ mb  v ■ i 
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0 0 10 
0 0 0 1 
0 0 0 0 
0 0 0 0 
0 0 0 0 
10  0 0 
0 10  0 


0 

0 


4.44 

-4.10 


6. OS  -1.03 


-0.73 

0 

0 


Figure  11.  Transformations  for  Plant  Canonical  Form  for  (A, 3,0 
in  Figure  9. 


For  the  reduced  order  observer: 


900  0 

_ 0 5625 

-33.3  0.66 

-1.96  -27.4 


R * diag [0 .0 1 , 0.01,  0.01,  0.01,  0.01] 


-1.14  1 0 
-0.05  0 1 


-37.9  36.8  524  0 0 

-3.39  -722  48. 3 0 0 


Other  design  parameters: 


1230  0 

0 900 


.515  0.00198  -0.069 

R = 0.00198  0.149  -0.104 

-0.069  -0.104  0.408 

V = 0 


Resulting  feedback  gains: 


0.08 

0.06 

2.23 

1.69 

-1.37 

2.14 

0.14 

4.57 

-3.23 

5.73 

-90.3 

-72.3 

0.74 

1.33 

-0.74 

2.89 

-0.56 

81.7 

10.2 

0.12 

2. 15 

Figure  13.  Design  Parameters  and  Resulting  Feedback  Gains  for  Sensitivitv 
Design  of  the  Aircraft  Control. 
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cas  possible  using  Che  solution  to  the  optimal  control  problem 


minimize  J = — u „x 1 (t)QK (t)  +u ' (t)Ru (t ) ]dt . 


The  arrangement  of  the  state  variables  was  as  in  (72).  It  is  well  known 
that  the  state  feedback  which  solves  (73)  is  given  by 

u(t)  m -R  ^B'Kx(t)  (7 

where  K is  the  nxn  positive  definite  solution  to  the  algebraic  Riccati 


equation 


KA  + ATK  + Q - KBR* Vk  = 0. 


The  matrices  Q,  R,  and  resulting  state  feedback  (74)  are  given  in 
Figure  14. 

As  is  conventionally  done,  the  state  feedback  matrix  was 
decomposed  into  the  five  columns  which  multiply  the  known  states  and  the 
two  columns  which  multiply  the  estimated  states.  The  resulting  feedback 
gains  used  in  the  system  with  the  reduced  order  observer  used  in  the 
sensitivity  design  is  given  in  Figure  15. 

As  an  evaluation  of  the  sensitivity  properties  of  these  two 
design  procedures.  Figures  17  to  20  show  the  trajectories  of  the  states  whie! 
were  to  converge  to  set  points  under  nominal  and  perturbed  conditions. 

The  set  points  and  perturbed  initial  conditions  are  given  in  Figure  lo. 

The  clearest  demonstration  of  the  improved  sensitivity  charac- 
teristics of  the  sensitivity  design  algorithm  is  a comparison  of  Figures  13 
and  20.  Clearly  the  Riccati  design  is  much  more  sensitive  to  the  Pertur- 


bations used.  The  solutions  presented  here  are  not  intended  to  represent 


A 


particularly  desirable  solutions  to  the  aircraft  control  problem,  but  to 


merely  demonstrate  the  sensitivity  design  procedure  under  some  practical 
design  conditions. 


i 


diagfl .44 , 

0.0025, 

12.25, 

12.25,  25 

, 0.0025, 

.01] 

diag [2 , 2, 

2] 

-2.31 

-0.02 

-0.14 

-2.47 

-0.10 

0.02 

-0.04 

0.02 

3.24 

2.46 

-0.13 

-0.40 

-1.32 

0.53 

-0.04 

-5.36 

-0.27 

0.09 

-3.51 

0.53 

-2.21 

Figure  14.  Design  Parameters  and  Resulting  State  Feedback  Gains,  F , 
for  Riccati  Design  of  the  Aircraft  Control.  R 
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-2.29 
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o 
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-0.35 

-2.47 

-1.04 
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0.53 

- 1 .46 

16.91 

22.98 

-0.13 

-0.40 

0.53 

-2.21 

0.47 

-63.53 
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0.09 

3.51 

Figure  15.  Feedback  Gains,  Used  to  Implement  the  Riccati  Design 
with  the  Reduced  Order  Observer. 
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Figure  16.  Set  Points,  E,  and  Initial  Condition  Perturbations, 
5x(0‘),  for  the  Aircraft  Control  Example. 
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6.  CONCLUSIONS 

The  advantageous  aspects  of  the  design  of  sensitivity  reducing 
dynamic  compensators  through  the  use  of  state  observers  have  been 
explained  in  the  body  of  this  thesis.  Some  of  the  shortcomings  of  the 
design  procedure  as  it  is  developed  here  should  be  pointed  out. 

Experience  bears  out  that  it  is  difficult  to  attain  desirable 
control  objectives  through  the  free  parameters  at  the  disposal  of  the 
designer.  This  is  because  of  the  complexities  of  the  equations  which  lead 
to  the  final  feedback  gains.  Although  the  observer  design  is  arbitrary 
and  determines  some  of  the  overall  feedback  system  eigenvalues,  the 
remaining  eigenvalues  are  not  easily  placed.  Consequently,  the  value  of 
the  algorithm  could  be  enhanced  by  further  research  of  the  dependence  of 
the  final  result  upon  the  free  parameters. 

More  insight  into  possible  applications  of  the  design  method 
could  be  gained  by  perhaps  considering  some  special  cases.  The  sensitivity 
weighting  matrix  H (65)  would  be  more  malleable  under  the  design  con- 
straints if  some  of  the  components  in  the  output  trajectory  variations  in 
(67)  and  (69)  were  not  present.  This  would  be  the  case  if  <*r.t=m.  , for 
instance,  which  would  often  occur  in  applications. 

In  conclusion  it  should  be  pointed  out  that  although  the  design 
procedure  appears  to  lead  to  a solution  of  the  regulator  problem,  it  has 
much  wider  application.  Any  design  which  calls  for  stabilising  compensa- 
tion may  have  improved  sensitivity  characteristics  by  employing  the  algorithm 
developed  in  this  thesis.  For  instance,  the  example  presented  in  the  last 
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APPENDIX 


DERIVATION  OF  THE  SENSITIVITY  REDUCTION 
CRITERION  FOR  THE  COUNTEREXAMPLE 


Consider  the  closed-loop  conf iguration  in  Figure  A.l  where  all 


variables  and  transfer  functions  are  scalars 


transfer  function  for  nominal  values  of  the  plant  parameters.  G(s)  and 


H(s)  are  the  controller  transfer  functions.  Let  6\’  (s)  and  5 Y (s)  be 


the  Laplace  transforms  of  the  first  order  output  variations  for  the  closed 


loop  control  and  the  nominally  equivalent  open-loop  control,  respectively 


resulting  from  differential  changes  in  the  plant  parameters 


It  is  shown  in  [4]  that 


where  S(s),  the  comparison  sensitivity  function,  is  defined  as 


Hence,  a necessary  and  sufficient  condition  for  a particular  closed-loop 


control  to  be  sensitivity  reducing  is  [4] 


The  example  of  Section  3.2  can  be  put  into  the  configuration  of 


Figure  A.l  by  defining  the  feedback  transfer  functions 


si 


where  M(s)  is  the  Laplace  transform  of  m(t)  defined  as 

A. 6.  m(t)  = 2.88xL(t)  - 10.18x?(t)  (A. 6) 

2^  and  2?  are  defined  as  the  observer  gains  in  equation  (20).  In  other 
words,  m(t)  is  the  feedback  signal  resulting  from  y(t)  = x2(t);  x (t), 
the  first  state  of  the  observer  (19);  and  the  optimal  feedback  gains  (17). 

Figure  A. 2 is  a block  diagram  of  the  example  in  terms  of  H^(s), 


Hn(s)  and  P^(s),  where 


P (s)  = -V 


s"-3s  - 2 


(A.  7) 


Figure  A. 3 is  an  equivalent  representation  of  the  feedback  system  in  a 
form  similar  to  Figure  A.l.  Hence,  for  the  example,  we  have  from  Figure  A. 3 
and  equation  (A. 2) 


S (s ) = 


Pn(s)  H:(s) 

1 -H^s) 


(A. 3) 


Substituting  A. 4,  A. 5 and  A. 7 into  A. 8 gives 


s'+-(8.38+2„>  s +(28.04-1.38  +5 .882  )s“- (25  .44-5 .642,  +6 .642,  )s 

i c 12 


-24. S+3. 762, -5.762, 

S(s)=- 3 2 . 

s +(1.30-*o)s  -(12.68-2  94. 302o)s^-(15. 26-4. 302,+2.222.,)s 
-4.44+2 . 222^ 

Equation  A. 9 in  conjunction  with  the  sensitivity  reduction  criterion  A. 3 
gives  the  inequality  (21)  in  section  3.2. 


(A. 9) 


